Digital Low-Pass Filter Implementation #1 



Define sampling frequency 
Fs :=4410C 
Define cutoff frequency 
Fc := 1400C 

Normalized cutoff frequency 

fc:=2^ 
Fs 

fc = 0.635 
Define Filter Q 
q :=10 

q 

Define Filter Parameters 



F:=2-sin|fc-— 

F= 1.68 
. . F 2 



G(x) 




arg(gain(D,x)) 



5000 




A =1.411 

B:=A 

Define Transfer Function (1 st column is numerator coefficients in transfer function 
and 2nd column is denominator coefficients in transfer function). So this filter is a 2nd 
order IIR. 



D:= 



D = 



B 1 

B -(2-FQ-F 2 ) 
0 (1-F-Q) 



1.411 1 
1.411 0:991 
0 0.832 



Plot frequency and phase response of the fitter 
x:=0,.00L .5 
G(x) :=| gain(D.x) | 



FIG. Al 




Di gital Hl-Pass Filter #1 



Define sampling frequency 
Fs :=4410C 
Define cutoff frequency 
Fc :=2000 

Normalized cutoff frequency 

fc .-2- — 
Fs 

Define Filter Q 
q:=2 

Q:=I 

q 

Define Piter Parameters 
F:=2-smlfc-— 

F = 0.284 

y 

2 

A = 0.04 



A :=- 



B:=A 



G(x) 1 




arg(gain(D,x)) 



5000 



'• 10 \44100 1 - 5 * 10 



2*10, 




Define Transfer Function (1st column is numerator c^fficients In f n f'^i on 2nd 
and 2nd column is denominator coefficients in transfer function). So this filter is a #ia 



order IIR. 
1 



D := 



-2 -(2-F-Q-F 2 ) 
1 (1-F-Q) 



D = 



1 1 

-2 -1.777 
1 0.858 



Calculate gain compensation 
4 

PK:=. 



(4- 2-F-Q-F 2 ) 

K:=-L 
PK 

K = 0.909 

Plot frequency and phase response of the filter 

x :=0, .001.. .5 

G(x) :=K-| gain(D.x) | 
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Digital Low-Pass Filter Implementation #2 (one-zero) 



Define sampling frequency 
Fs :=4410C 
Define cutoff frequency 
Fc :=2000 

Normalized cutoff frequency 

fc:=2.^ 
Fs 

fc= 0.091 
Define Filter Q 
q :=10 

Q.--L 

q 

Define Filter Parameters 
n 



G(x) 




arg(gain(D,x)) 



5000 



1-10* ,., m 1.5«10 4 
x -44 100 



F:=2-sin^fc 
F = 0.284 

2 

A = 0.04 
B:=A 

Define Transfer Function (1st column is numerator coefficients in transfer function 
and 2nd column is denominator coefficients In transfer function). So this filter is a 2nd 
order IIR. 




D:= 



D = 



B -(2-F-Q-F 2 ) 
0 (1-F-Q) 

0 1 
0.04 -1:891 
0 0.972 



Plot frequency and phase response of the filter 
x:=0,.00L. .5 
G(x) :=| gain(D.x) | 



FIG. A3 



Digital Shelving EQ section based on an All-pass filter section 



Define sampling frequency 

Fs :=4410C 

Define critical frequency 

Fc := 1000C 

Normalized critical frequency 

„ Fc 
wc :=2-it • — 
Fs 

Define gain/cut 
k:=4 

Define Filter parameters 



G(x) 



gamma. = 



tanfrl-1 



tan| — | -t- 1 



\ 2 



o.i 



Define Transfer Function (1 st column is numerator coefficients in transfer function 
and 2nd column is denominator coefficients in transfer function). So this filter is a 1 st 
order IIR. 



A := 



-gamma 1 
- 1 gamma 

0.073 1 
-1 -0.073 



Plot frequency and phase response of the filter 

x:=0,.00L0.5 

G(x) :=| gain(A.x) | 



arg(gain(A,x)) 



To produce desired parametric EQ embed all-pass section in below structure 

1+K 

, >C X 




in a DSP implementation, the filter output is computed as per the above signal 
flow. Below we calculate the overall transfer function of the above system by 
transforming the numerator of the all-pass section as necessary. 



FIG. A4 (Sheet 1 of 2) 



D:=A 



G:=D 1 extract denominator 



c ._r> <0> 

r.-U extract numerator 




FIG. A4 (Sheet 2 of 2) 0006 




Digital Notch/Parametric EQ section based on an all-pass section 



Define sampling frequency 
Fs :=4410C 

Define critical frequency 
Fc:=5500 

Normalized critical frequency 

wc :=2-Jt • — 
Fs 

Define Filter Q 
Q:=4 

Define gain/cut 
k:=4 

Define Filter parameters 
gamma :=-cos(wc) 

wc 



G(x) 



fbeta :=- 

2Q 

fbeta :=ififbeta>— ,— , fbeta 
\ 4 4 

,_ I- tan(fbeta) 

DCul • " 

1+- tan (fbeta) 

Define Transfer Function (1st column is numerator c 
and 2nd column is denominator coefficients in transl 
order IIR. 

beta 1 
A := gamma-(l-t-beta) gamma(l-f-beta) 
1 beta 

0.821 1 
A= -1.29 -1.29 
_ 1 0.821 

Plot frequency and phase response of the filter 
x:=0,.OOL 0.5 
G(x) :=| gain(A.x) | 



arg(gain(A.x)) 0 



"2 



~4 



0.1 



0.2 



0.1 



0.2 



To produce desired parametric EQ the all-pass section is embedded in the below structure: 

1+K 




In a DSP implementation, the filter output is computed as per the above signal 
flow. Below we calculate the overall transfer function of the above system by 



FIG. A5 (Sheet 1 of 2) 




transforming the numerator of the all-pass section as necessary. 

D:=A 

r—n <l> 

<J •- u extract denominator 

n-n <0> 

r.-D extract numerator 




Magnitude Response: 



FIG. A5 (Sheet 2 of 2) 
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Digital Implementation of a Weighted Notch #1 

Consists of 3 parametric EQ section designed using the procedure described in 
"Notch/Parametric EQ section based on an all-pass section" 

2 low Q gain section surround the critical frequency notch. 

For compensation section #1 and #2 (D1 and D2 respectively) Q is 2 and k is 1 .5 

Stag in the below transfer functions. The resulting filters are two 2nd order IIP, filters where column 0 represents the 
numerator and column 1 represents the denominator. 



Dl := 



D2:= 



1.076 1 
- 1.283 -1.283 
0.619 0.695 

1.089 1 
-1.079 -1.079 
0.554 0.643 



For critical frequency section #3 (D3 below) Q is 10 and k is 0. The resulting filter is a 2nd order IIP. filters where column 0 
represents the numerator and column 1 represents the denominator. 

"0.962 1 
D3:= -1.363 -1.363 
0.962 0.925 , 

Plot of the overall frequency response of the 3 IIP. filters cascaded 

X :=0,.001..0.5 
Gl(x) :=| gain(Dl.x) | 
G2(x) :=i gain(D2,x) | 
G3(x) :=| gain(D3,x) | 
H(x) :=201og(Gl(x) G2(x) G3(x)) 



H(x) 




-10 - 



FIG. A6 
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Digital Implementation of a Weighted Notch #2 

Consists of 3 parametric EQ section designed using the procedure described in 
"Notch/Parametric EQ section based on an all-pass section" 

One medium Q gain compensation section centered around two notches. 

For compensation section #1 (D1 below) Q is 4 and k is 4 resulting in the below transfer function. The resulting filters is a 
2nd order MR filters where column 0 represents the numerator and column 1 represents the denominator. 



Dl:= 



1.268 1 
-1.29 -1.29 



0.553 0.821 

For critical frequency notches sections #2 and #3 (D2 and D3 below respectively) Q is 10 and k is 0. The resulting filters 
are two 2nd order IIP. filters where column 0 represents the numerator and column 1 represents the denominator. 



D2:= 



0.963 1 
-1.383 -1.383 
0.963 0.926 



0.962 1 
D3:= -1.343 -1.343 
0.962 0.923, 

Plot of the overall frequency response of the 3 MR filters cascaded: 



x:=0,.OO0L.0.5 

Gl(x) :=| gain(Dl.x) | 

G2(x) :=| gain(D2,x) j 

G3(x) :=| gain(D3,x) | 

H(x) :=201og(GKx) G2(x) G3(x)) 



(cascaded response) 



H(x) -20 - 




-40 - 



-60 



FIG. A7 
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FIG. 8a 
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FIG. 8c 
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FIG. 8d 
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FIG. 16 
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FIG. 3a 
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FIG. 3b 
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FIG. 3c 
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